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Abstract 

The damping a of the spinwave resonances in 75 nm, 120 nm, and 200 nm -thick Permalloy 
films is measured via vector-network-analyzer ferromagnetic-resonance (VNA-FMR) in the out-of- 
plane geometry. Inductive coupling between the sample and the waveguide leads to an additional 
radiative damping term. The radiative contribution to the over-all damping is determined by 
measuring perpendicular standing spin waves (PSSWs) in the Permalloy films, and the results are 
compared to a simple analytical model. The damping of the PSSWs can be fully explained by 
three contributions to the damping: The intrinsic damping, the eddy-current damping, and the 
radiative damping. No other contributions were observed. Furthermore, a method to determine 
the radiative damping in FMR measurements with a single resonance is suggested. 


Contribution of NIST, not subject to copyright 


1 



I. INTRODUCTION 


Excited magnetic moments relax towards their equilibrium orientation due to damp¬ 
ing. Several physical mechanism can cause damping. Many mechanisms, like eddy current 
damping in conducting ferromagnets, were already identihed in the 1950s. More recently, 
enhanced damping due to spin pumping from a ferromagnet into an adjacent metallic layer 
was identihed, and remains a topic of ongoing invest igatiorP®. Furthermore, wavenumber- 
dependent contributions to the damping caused by intralayer spin pumping have been theo¬ 
retically predictecP® and currently are the subject of experimental invest igatioiP^. Another 
damping process, referred to as radiative dampin^^I^, has been known to exist since the 
1970s and is purely due to inductive coupling between the sample and the waveguide in ferro¬ 
magnetic resonance (FMR) experiments. More recently, this phenomenon has been further 
investigated in the context of strong magnon-photon coupling experiments, with possible 
applications in quantum information processin^^^. In these quantum-coherent experiments, 
radiative damping was identihed as a manifestation of non-resonant magnon-photon cou¬ 
pling, and it was determined that such couplings! jg indeed a source of extrinsic line width 
in cavity-based FMR studies. In a radiative damping process, the time-varying magnetic 
hux associated with the dynamic magnetization generates microwave-frequency currents in 
the proximate conductor of a microwave waveguide that carries the resultant power away 
from the sample. This process is similar to that exploited in eddy-current brakes and can 
be seen as a non-local counterpart to the eddy-current damping in conductive ferromagnets. 
To determine the magnitude of radiative damping in magnetic thin hlms, we used broad¬ 
band vector-network-analyzer ferromagnetic resonance (VNA-FMR) to measure damping in 
Nio. 8 Feo .2 Permalloy (Py) hlms with thicknesses 6 varying between 70 nm and 200 nm. By 
use of the geometry sketched in Fig. [^a), we determine the total damping for each mode 
On as a sum of intrinsic damping a™*, eddy-current damping and radiative damping 
We then perform a quantitative analysis of the PSSW resonance helds, amplitudes and 
damping to extract the diherent contributions to We hnd that eddy current damping is 
only signihcant for the lowest order mode, the radiative damping strongly ahects the hrst 
hve modes, and no additional contributions to the damping are detectable for spin waves up 
to k = 1.75 X 10®cm“^. This last hnding is in contrast with reports of exchange mediated 
damping in both nanostructure^ and thin hlm^. 
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II. DAMPING MODELS 


According to Faraday’s Law, the time-varying flux of a processing magnetic moment 
generates an ac voltage in any conducting material that passes through the flux. As shown 
in Fig. [TJ spin wave precession in a conducting ferromagnet on top of a coplanar waveguide 
(CPW) induces ac currents both in the ferromagnet and the CPW. The dissipation of these 
eddy currents in the sample, and the flow of energy away in the CPW give rise to two 
contributions to magnetic damping. Historically, the damping caused by eddy currents in 
the ferromagnet is called eddy current damping, while the damping caused by the eddy 
currents in the waveguide is called radiative damping 

Eddy current damping has been recognized since the 1950^^^. For the lowest order mode 
in FMR, 

( 1 ) 

16 p 

with the resistivity p, saturation magnetization Mg, the vacuum permeability po, the gyro- 
magnetic ratio 7 , and the sample thickness 5 (see derivation in Appendix Sec. B). We 
introduce a correction factor C to account for details of the eddy current spatial profile. As 
shown in a later section, for all higher order PSSW modes investigated in this study 

is much smaller than that of the lowest order mode. 

We now turn to the radiative damping We consider the experimental geometry 

sketched in Fig. (a). A ferromagnetic sample with thickness 6 and length I is placed on 
top of the center conductor of a coplanar waveguide with width W. The sample dimension 
along X is much larger than W. The sample and CPW are separated by a gap of height d. 
An external dc magnetic held Hq is applied perpendicular to the sample plane, and the spin 
wave resonances (SWR) are driven by microwaves in the CPW at resonance frequency /. A 
fraction of the ac magnetic induction B due to the dynamic component of the magnetization 
I] z) wraps around the center conductor. 

To derive a quantitative expression for we start by calculating the x 

component of the driving held that is generated by an excitation current I in the 
center conductor. We assume M{r, x, z) is uniform along y, but we allow for variation along 
X and 2 . To estimate h^{I-,x,z) we use the Karlqvist equatiorP^ 


h^{I;x,z) = 


27iW 


arctan 


x + W/2\ 


arctan 
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W/2 
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Figure 1: Schematic of the radiative damping process, (a) M is the dynamic magnetization, Hq 
the applied external held and B is the magnetic inductance due to the x-component rux of the 
dynamic magnetization. W is the width of the center conductor, I the length of the sample on the 
waveguide, 5 the thickness of the sample and d the spacing between sample and wave guide, (b) 
Simplihed depiction of the PSSW eigenfunctions Qn for mode numbers n= 0, 1, 2, 3. We exemplarily 
used boundary conditions that are completely pinned on one side and completely un-pinned on the 
other side. The origin of the coordinate system is indicated. 
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This microwave field can excite PSSWs in the sample. Schematic mode profiles for the 
fnndamental mode {n = 0) and the first three PSSW modes are shown in Fig. [^b), where 
we use unpinned boundary conditions at the top surface and pinned boundary conditions at 
the bottom surface. As shown in Fig. [^b), the mode profiles describe a z-dependence of the 
dynamic magnetization components and m^. In the perpendicular geometry used here, 
|m*| = \m^\ everywhere, i.e. the precession is circular. In what follows, we will only discuss 
m*, the dynamics of which are inductively detected in the measurement. For a PSSW with 
mode number n, fh^{x,z) = qn{z)xn {<ln{z)h^{r,x, z)) where () denotes spatial averaging 
in X and 2 ; directions, as defined in the Appendix, Xn = is the diagonal component of 
the magnetic susceptibility of the n-th order mode, and —1 < qn{z) < 1 is the normalized 
mode profile (eigenmode), an example of which is sketched in Fig. (b) for n = 3. The 
mode inductance L„ is given by = XnLn, where, as detailed in the Appendix, we define 
a normalized mode inductance for the nth PSSW mode. 


Ln = ^{qn{z)h^{I;x,z)fW5. (3) 

Ln, as explained in the appendix, no longer has any dependence on magnetic field or 
excitation frequency. In the simplest case of a uniform magnetization profile qo{z) = 1 
(FMR-mode) and uniform excitation field h^{r,x, z) = h*(/; 0, 0) = I/{2W), the normalized 
inductance is Ln = /io^//(4IF). 

The x-component of the dynamic magnetization rrin^x, z) produces a net flux = Xn^Ln 
that threads around CPW center conductor, leading to a power dissipation 

= (4) 

where Zq is the waveguide impedance (in our case Zq = 50 hi) and uj is the angular frequency 
of the magnetization precession. With Eq. (|^, the power dissipation rate = Pn/En 

can be calculated, where En is the energy of the dynamic component of the magnetization 
derived in the Appendix. This power flow from the sample to the waveguide leads to the 
radiative damping contribution 
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where rj = 5/ {a dz\qn{z)f'^ is a dimensionless parameter that accounts for the actual 
mode prohle in the sample, see Appendix Sec. A. In the case of sinusoidal PSSWs, rj = 1/2, 
and for a completely uniform mode prohle, i.e. qn{z) = 1, rj = 1/4. From Eq. ([^, it is 
evident that is proportional to for u > 0. In the simplest case of uniform driving 
held = I/{2W), the radiative contribution is given by 

° Zo 2ZohF ^ ^ 

Note that the radiative damping thus depends on the sample and waveguide dimensions, in 
particular linearly on the sample thickness. Unlike eddy-current damping, is indepen¬ 
dent of the conductivity of the ferromagnet, hence this damping mechanism is also operative 
in ferromagnetic insulators. 


III. SAMPLES AND METHOD 

We deposit Ta(3)/Py((5)/Si3N4(3), Ta(3)/Py(h)/Ta(5), Ta(3)/Py((5) and Py(5) layers on 
100 pm thick glass substrates by DC magnetron sputtering at a Ar pressure of 0.7Pa (~ 
5 ■ 10“^Torr) in a chamber with a base-pressure of less than 5 ■ 10“® Pa (~ 4-10“®Torr); 
where 6= 75 nm, 120 nm and, 200 nm is the Permalloy thickness. The Py thickness was 
calibrated by x-ray rehectivity. We estimate that the damping enhancement due to spin 
pumping into the Ta layer is two orders of magnitude smaller than the intrinsic damping of 
the Permalloy layer for Permalloy samples of these thicknesses. The various combinations of 
capping and seed layers are chosen to determine the sensitivity of our results on the spinwave 
boundary conditions and the resultant mode prohles. Prior to deposition, the substrates are 
cleaned by Ar plasma sputtering. The samples are coated with approximately 150 nm of 
PMMA in order to avoid electrical shorting when samples are placed directly on the CPW. 
The CPW has a center conductor width of hF = 100 pm. The SWR are characterized using 
held-swept VNA-FMR?^® in the out-of-plane geometry (see Fig. with an external static 
magnetic held Hq applied perpendicular to the sample plane. The excitation microwave held 
h^{x, y) is applied over a frequency range of 10 GHz to 30 GHz. A VNA is used to measure 
the complex 521 transmission parameter (ratio of voltage applied at one end of the CPW 
to voltage measured at the other end) for the waveguide/sample combination. The change 
in ^21 due to the FMR of the sample is then htted with a linear superposition of complex 
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susceptibility tensor components Xn, 


N 

AS^^{Hq) = AnXn{Ho)e^'^" + linear background (7) 

n=0 


with the mode number n, phase (j)n, and dimensionless mode amplitude An, as dehned in the 
Appendix. A complex linear background and offset is included in the £t. The susceptibility 
components are derived from the Landau-Lifshitz equation for the perpendicular geometry; 
in the fixed-frequency, swept-field configuration, we obtairP^ 


Xi^o)n 


_ M,{Ho - Mf - _ 

{Ho - Mf - HzY - - lAHniHo - Mf - if-) 


( 8 ) 


with if®® = uj/ ( 7 /io) and = Mg — Hk, where ffk is the perpendicular anisotropy, is 
the exchange field (dehned below), and Aff„ is the linewidth. An example of the resulting 
fits for the complex S '21 data is shown in Fig. (a) and (b). 


IV. EXPERIMENT 

We detect both even and odd PSSW modes. If we assume a uniform excitation held 
and Dirichlet boundary conditions (completely pinned), only odd modes would be detected. 
Alternatively if we assume Neumann boundary conditions (completely unpinned), only the 
fundamental mode would be detected. 

Two ehects can contribute to our ability to detect all the PSSW modes. First, the 
excitation held prohle might not be uniform due to eddy current shieldin ^^^P^l Second, the 
interfacial boundary conditions might be asymmetrical, as alluded to above. According to 
the criterion in Ref. y , the threshold sheet resistance for the onset of eddy current shielding 
at 20 GHz is 0.065 G/D. We estimate that the sheet resistance for our 200 nm is in excess of 
0.345 G/n, so we conclude that the eddy current shielding is relatively weak for our samples. 

On the other hand, all modes are in principle detectable if we assume asymmetric in¬ 
terfacial anisotropy. For the sake of simplicity of the analysis, we will assume interfa¬ 
cial anisotropy for a single interface, then use an optimization approach to determine the 
wavenumber of the modes that is consistent with such a hypothesis. However we must 
emphasize that this approach does not provide a unique ht for the measured distribution 
of resonance helds for the PSSW spectrum, but simply allows us to accommodate for the 
wavenumber values required to be consistent with the measured spectrum. As such, the 
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Figure 2: Measured ^21 transmission parameter (black circles) at 20 GHz and the mnlti-peak 
-susceptibility fit (red line) for the (a) real part and (b) imaginary part obtained with the Ta(3)- 
Py(200)-Si3N4(3) sample. The hrst 6 modes are shown, (c) The exchange held (black squares) 
and exchange held ht, from Eqs. (§ and (red crosses) for all 13 detected modes plotted as a 
function of the htted wave numbers kr,. 


fitted value for iTg is to be interpreted as no more than a self-consistent value associated 
with only one of many possible scenarios. 

If we assume negligible magnetocrystalline perpendicular anisotropy hfk, is related 
to the exchange held via 


i/r = +Ms, 

2 

ki. 


withi7f= 

hoM 


(9) 


Here, kn is the spinwave wavevector, and Hex is the exchange energy that is related to the 
spinwave stiffness D via D = Qj^ other hand, if we want to include interfacial 

anisotropy for a single interface in our analysis, we can numerically solve the transcendental 
equationl^ 


--kna+ + 1 ) 

2 2A,^kn ) 
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where is the interfacial anisotropy, and a = 0.3547 nm is the lattice constanlP^. We 
minimize the residue of the £t of Eq. ([^ to with the fitting parameters M^, Aex, and 
Kg from Eq. (10) by use of a Levenberg-Marquardt optimization algorithm. This yields the 
pairs {kn,H^) shown in Figj^ (c) for all modes. 

From the fit, we obtain a saturation magnetization of = 1.02 ± 0.01 T, in agree¬ 

ment with that determined by magnetometry. The exchange stiffness constant oi D = 
3.22 ± 0.04meVnm^ is close to a value of H ~ 3.1 meVnm^ reported by Maeda et al^. 

The exchange £t also yields a single surface anisotropy that depends on the 
cap and seed layer conhgurations. For the Ta(3)-Py((5)-Si3N4(3) sample series. 
Kg, = (5.1 ±0.8) X 10“^J/m^, while all the other samples have a higher of 
(7 ± 1) X 10“^J/m^. All values for Ks are in the range of other reported interface 
anisotropies for Permalloy layers of these thicknesse^. 


We now turn to the linewidth Ai7„ and the amplitude for the individual modes. The 
Gilbert damping parameter is extracted from the slope of the linewidth vs. frequency / 
plol!^ shown in Fig. [^a) via 

± ( 11 ) 

where Ai7° is the inhomogeneous broadening that gives rise to a nonzero linewidth in the 
limit of zero frequency excitation. The normalized inductance of the modes is extracted 
in a similar fashion from the dependence of the mode amplitude on the frequency /, see 
Fig. I^b) and Eq. (38) in the Appendix 


A„ = 27r/^±A°, (12) 

where A\ is an offset for each mode. A° is a phenomenological fitting parameter, which is 
not yet fully understood. 

We plot an and as a function of mode number n in Fig. (c). The damping and 
the normalized mode inductance are found to be proportional. In order to explore this 
correlation, we plot vs. in Fig. |^a). Here, the data for vs. are linearly 
correlated for all modes except for n = 0, as seen by the linear £t (line) to the data for 
77, > 1. This is as expected for the radiative damping model, as summarized in Eq. ([^. The 
additional damping of the fundamental mode is interpreted as the result of eddy current 
damping, as quantihed in Eq. Q. In Fig. igb), we plot the residual Aa„ of the linear £t 
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Figure 3: Parameter extraction for the first 5 PSSWs of the Ta(3)-Py(200)-Si3N4(3) sample, (a) 
Extraction of a from the linewidth (data points) via linear fits (lines); staggered for display, 

(b) Extraction of the normalized mode inductance (data points) from the resonance amplitude 
An{f) via linear fits (lines), (c) an (black squares) and (red crosses) for each PSSW. 


shown in Fig. |^a) for all modes. Aa„ is negligible for all modes except for n = 0. We 
extract Aa„=o for all the samples and plot Aan=o vs. <5^ in Fig. |^c). 

It appears that Aq!„=o for all the samples scales linearly with 6^, as expected from Eq. ([^ 
for eddy current damping. Simultaneous weighted hts of all the data to Eq. ([^ yields 
C = 0.4 ± 0.1. This value suggests a localization of eddy currents, since C corrects for the 
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Figure 4: Damping an and inductance Ln for the Ta(3)-Py(200)-Si3N4(3) sample, (a) Linear fit 
of a to Eq. ([^ where the fit is constrained to the n= 1, 2, 3, 4 modes, (b) The residual of the 
linear fit, showing enhanced damping for the 0-th order mode (black). We attribute the enhanced 
damping to an eddy current contribution, (c) Enhanced 0-th order mode damping for all samples. 
The red line is a fit of the data points to the eddy current damping model from Eq. Q. 


eddy current distribution in the sample. 

For the n > 1 modes, it can be showrpS that oc 1/A;^. The calculated wavevectors 


from Eq. (10) for the n = 1 mode of all the samples is at least a factor three larger than that 
of the n = 0 mode and, therefore, the eddy current damping of the n = 1 mode is predicted 
to be approximately one order of magnitude smaller than the eddy current damping of the 
n = 0 mode. Thus, the eddy current damping of the n > 1 modes is negligible to within 
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the error bars, i.e., ~ 0 for n > 1. This supports the analysis of the data in Ref. [gj], 

which also neglects the eddy current damping in higher order modes. 


V. EXTRACTION OF THE RADIATIVE CONTRIBUTION TO THE DAMPING 


By use of Eq. ([^ and our htted value of (7 = 0.4, we subtract the eddy current contri¬ 
bution to the damping of all the u = 0 modes to obtain a corrected damping value 0^=01 
where = otn=Q — = 0.4). The corrected data for all the modes are plotted in 

Fig.|^ 

Figures (a) to (c) group all data obtained for a set of samples with identical Py thickness 
5. The lines are linear hts to Eq. ([^. For each thickness 5, we observe a signihcant correlation 
of On and for all seed and cap layer conhgurations, as expected for a radiative damping 
mechanism. 

Furthermore, by use of Eq. (§ for the n = 0 mode of the 75 nm thick sample, using a 
value oi 7 ] K. 0.46 as determined in the Appendix, we estimate ~ 0.00023 

The experimentally determined value is ~ 0.00035 ± 0.0001. The deviance from 
the calculated value is possibly due to non-uniformities of both the excitation held and 
magnetization prohle in Eq. (§, that requires the solution of the integral in Eq. ( [2^ . 
Nevertheless the estimated value for is of the correct order of magnitude. 


We determine the intrinsic damping from the = 0 intercept of the linear hts in 
Fig. g We plot a™* for the three values of 6 in Fig. (right scale). We hnd that a™* is 
approximately constant to within ±5% for all samples. In addition the average value over 
all the him thicknesses is in reasonable agreement to the previously reported value of 
0.006 (dotted red line)l^. 

The other htting parameter rj, extracted from the slope of vs. Ln, is also plotted in 
Fig. 1^ (left scale). For anti-symmetric boundary conditions, rj = 1/2 is expected, whereas 
for the uniform mode, p = 1/4. 

We see that the htted values lie exclusively within these extremes, within error bars. 


There have been recent reports of a non-zero, wavenumber-dependent component for 
damping for both localized eigenmodes in magnetic nanostructure^ and PSSWs in thick 
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Figure 5: Dependence of damping a on the normalized mode inductance after correction for 
the eddy current damping Aoq of the fundamental mode, plotted for all sample configurations 
and different thicknesses: (a) 5 = 200 nm, (b) 5 = 120 nm and (c) 5 = 75 nm. The red lines are 
weighted linear fits to the data by use of Eq. Q, that describes the radiative component of the 
damping. 
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Figure 6: Mode profile parameter rj (black squares, left axis) and intrinsic damping a\^t (red circles, 
right axis) as a function of Py thickness 5. The mode profile parameter rj lies between the value 
of 1/2 for sinusoidal PSSWs with anti-symmetric boundary conditions (dashed black line) and the 
value of 1/4 for the uniform mode (both dashed black lines). The intrinsic damping is close to 
CKint = 0.006 (dotted red line). 


Permalloy film^. Such exchange-mediated damping of the form := was originally 

predicted by Baryakhtar based on symmetry alon^. Nembach, et al.^, obtained a value of 
^ex = 1.4 nm“^, whereas Li, et al.0 found a much smaller value of 0.09 nm To determine 
whether wavenumber-dependent damping is apparent in our data, we examined the residual 
damping after subtraction of both the intrinsic damping a™* and the radiative damping 
from all the modes, as well as subtraction of the eddy current damping from the n = 0 mode. 
The residual damping a’’®® is plotted in Fig. [^(b). Within the scatter of = ±0.001, W®® does 
not have any clear dependence on k. Thus, we obtain an upper bound of < 0.045 nm“^ 
for this particular system, given the sensitivity of our measurements. For comparison and to 
ensure that the subtraction of a™*, and q;®'^^^ did not hide a potential contribution 
the measured damping of the Ta(3)-Py(200)-Si3N4(3) sample up to the n = 10 mode is 
shown in Fig. (a). For n > 5 the measured damping scatters around the for the 200nm 
samples determined intrinsic damping a™* and no trend for higher mode numbers (larger k 
values) is discernible. 
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Tserkovnyak, et al., calculated the damping coefficient A^x in terms of a microscopic 
model for the diffusive transport of dissipative transverse spin current within a ferromagnetic 


metaP. The theory in Ref. y framed the exchange-mediated damping in terms of a so-called 
transverse spin conductivity ct_l, 



(13) 


where 



(14) 


with the exchange splitting huex, the conductivity a, the spin scattering time r, and 
transverse spin scattering time tj_. Given that hcUex ~ leV for Permalloy, the maximum 
value for Rgx predicted by the transverse spin current theory is 0.001 nm“^. Insofar as 
we are not able to observe any such wavenumber-dependent damping down to the level of 
0.045 nm“^, our results are consistent with the predictions of the microscopic theory. 

While the theory in Ref. [gj] is specihc to the microscopic mechanism of transverse spin 
accumulation in a metallic ferromagnet, the phenomenology of exchange mediated damping, 
as described in Ref. [g], is not limited to such a microscopic mechanism. As such, it remains 
plausible that extrinsic material-specihc parameters that have not yet been identihed could 
be responsible for the previously reported values for damping. For example the presence 
of anti-symmetric exchange at interfaces, i.e. the Dzyaloshinskii-Moriya interaction (DM1) 
could enhance the coupling between magnons and Stoner-excitations insofar as the DM1 
gives rise to exotic spin textured with nanometer length scales, that are comparable to 
the wavelength of low energy Stoner-excitationd^. Thus, the results of Ref. [loj] could be 
a manifestation of interfacial enhancement for A“, insofar as the magnetic dims used in 
Ref. [loj] are only 10 nm thick. 

In another experiment, we further validate the presence of radiative damping and demon¬ 
strate an alternative method to determine by varying the distance d in Eq. (|^ between 
the sample and waveguide. To this end, we insert a. d = 200 pm glass spacer between the 
sample and waveguide. By comparing h(0,0) to h(0,200 pm) via Eq. ([^, we estimate that 
the insertion of the spacer decreases the microwave magnetic held by about a factor of 6.25. 
Referring to Eq. ([^, the normalized mode inductance Ln decreases by a factor of = 40. To 
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Figure 7: (a) The measured damping for the first 11 PSSW modes of the Ta(3)-Py(200)-Si3N4(3) 
sample. The enhanced damping due to inductive coupling to the waveguide and eddy currents 
in the sample only affects the first five modes at wavevectors < 7 x 10®cm“^. (b) The residual 
damping for all detected modes for all samples is plotted against their respective wave vector k. 
Within the scatter, no dependence of the residual damping on k is observed. 

determine the effect of the reduced inductive coupling on the radiative damping, we used 
VNA-FMR to measure the first 4 modes for the Ta(3)-Py(120) sample with and without 
the spacer. The effect of the spacer can be seen in the raw data, reducing the linewidth 
of the hrst two modes measured at 10 GHz in the 120 nm samples by approximately 6 Oe, 
well outside error bars. The fitted values of are shown in Fig. (a). Indeed, de¬ 
creases on average for all modes by a factor of = 50 after inserting the spacer, in good 
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agreement with the predictions of Eq. ([^ and (|^. Thus, we will assume that is neg¬ 
ligible when the spacer is used. The data for the damping of the hrst four modes, both 
with and without the spacer, are plotted in Fig. (b). Indeed, the damping determined 
from the measurement with the spacer layer (circles) is consistently lower than that found 
without the spacer layer (squares). The line in Fig. (b) is the previously determined in¬ 
trinsic damping. Under the assumption that the radiative damping contribution is given by 
= an{d = 0) — an{d = 200 pm), we plot vs Ln{d = 0) in Fig. (c). The line is 
the calculated where we used Eq. ([^ with rj = 0.35 and 5 = 120 nm, as determined 
from the hts in Fig. Good agreement between the calculated and measured values for 
are obtained, which demonstrates the self-consistency of our analysis. Of great importance 
is that the spacer-layer approach can also be used to determine the radiative contribution 
to the damping in the absence of PSSWs (single resonance). By measuring a for varying 
distance d between sample and waveguide and extrapolating a to d —?• oo, both the intrinsic 
value for the damping and the radiative contribution can be determined, under conditions 
where eddy current damping is negligible. 


VI. SUMMARY 

In summary, we identihed three contributions to the damping in PSSWs: Intrinsic damp¬ 
ing a™*, eddy current damping and radiative damping The latter exhibits a 

linear dependence on the normalized sample inductance Ln in a waveguide based FMR 
measurement. We attribute this linear dependence to radiative losses that stem from the 
inductive coupling between the sample and the waveguide. The radiative damping term is 
inherent to the measurement process and is thus present in all FMR measurements. The 
radiative damping constitutes up to 40 % of the total damping of the spin wave modes 
in our 200 nm thick Permalloy Elms. Furthermore, the radiative damping can be already 
important for much lower him thicknesses, in materials with small intrinsic damping. 

As an example, the radiative damping calculated from Eq. (§ for a 20 nm thick and 
1 cm long sample of Yttrium-Iron-Garnet (YIG), measured on a 100 pm wide wave guide is 
^ 1.26 ■ 10“^ . When compared to the reported value for the damping of a = 2.3 ■ 10“^,^ 
we see, that the radiative part of the damping, among other^, can substantially inhuence 
the determination of amt- As such, careful analysis of a vs. inductance is required to isolate 
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Figure 8: Measurement of the first four PSSWs of the Ta(3)-Py(120) sample with and without a 
spacer inserted between sample and CPW. (a) Inductance determined for the sample directly 
on the CPW (black squares) and for a 200 p.m spacer between sample and CPW (red circles), (b) 
The resulting damping constants for both measurements (same symbols and colors). The red line is 
the previously extracted intrinsic damping (c) The difference between the damping with and 
without the spacer (black squares) is in good agreement with the radiative damping from Fig. [^b) 
(gray line). 


the radiative damping contribution. 
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VIII. APPENDIX 


A. Derivation of 


In this section, we derive model equations for the normalized mode inductance mode 
amplitude An, and the radiative damping We are restricting our analysis to the case 
of ideal perpendicular standing spin wave modes that only vary through the him thickness 
without any lateral variation. It is assumed that the excitations are in the perpendicular 
geometry with the magnetization saturated out of the him plane. As such, the response to 
the z-coordinate component of the microwave excitation held above the waveguide can be 
neglected. In addition, the magnetization precession is always circular. As such, the magne¬ 
tization dynamics in the x- and y-coordinates in response to the microwave held generated 
by the waveguide are degenerate, outside of a phase factor of 7r/2. This eliminates the need 
to explicitly consider the full Polder susceptibility tensor in the calculation of the sample 
response to the excitation held. The sample dimensions are I along the waveguide direction, 
6 in thickness, but inhnite in the lateral direction. 

We begin by introducing the concepts of a spin wave mode susceptibility Xn, and the di¬ 
mensionless, normalized spin wave amplitude gni^) for the nth spin wave mode, such that 
the magnetic excitation of amplitude in a:-direction that results from the ap¬ 

plication of a microwave magnetic held of amplitude h^{I; x, z), given in Eq. ([^, driven by 
an ac current I = Vin/Z q in an applied held Hq is given by 

fhn {Ho, I-,z) := ml {z) = qn {z) Xn {Hq) {qn {z) (/; x, z)) (15) 


where the quantity in brackets is simply the overlap integral of the excitation held and 
the spatial prohle of the nth spin wave mode. The magnetic excitation of amplitude in 
y-direction {z) can be written in a similar way. In the trivial case of a uniform excitation 
held and uniform spin wave mode, we recover the usual relation between the excitation 
held and the magnetization dynamics via the Polder susceptibility tensor component, x^^- 
However, if the product of the mode prohle and excitation held has odd spatial symmetry, 
dynamics are not excited, as we expect. The overlap integral is nothing more than the 
spatial average of the mode/excitation product: 


{qn{z)h^{r,x,z)) 


1 

W6 





dzqn{z)N^{I-,x,z) 


(16) 
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First, the power transferred to the waveguide via inductive coupling with the spin wave 
dynamics is given by 


where 


Pn 


mn (go,/)| 

2Zq 


2 


(17) 


oo 5-\-d 

dt^n {Ho, I) = fioi j dx j dz {dtml {z)) {x, z ), (18) 

—oo d 

with If (x, z) = h'^{P,x, z)/I. 

It is important to recognize at this point that the power dissipation is not constant with 
time, given that is proportional only to dtm^. As such, the damping associated with the 
re-radiation of the microwave energy back into the waveguide is best characterized with an 
anisotropic damping tensor, to be elaborated upon more fully later in this Appendix. To 
calculate the energy of the spin wave mode, we start by dehning a spatially averaged spin 
wave excitation densit}^. 


IZo If * (*)) (™i; («))* - (St™?, (*)) (™? (*))' 

4uiSW 




(19) 


We can then calculate the magnon density Nn associated with the nth spin wave excitation 
as 


= 




2giiBMs 

The total energy associated with the spin wave mode is given by 


iM, 

The energy dissipation rate (1/Ti)^ for the nth mode is therefore 


( 20 ) 


( 21 ) 


1 


Pn 

Er, 



oo S-\-d 

2 

2fiQiuJM 

f dx f dz{dtml{z))h^ {x,z) 


—oo d 



OO S-\-d 


, ( 22 ) 


f dx f dz (z)) {ml {z})* - {dtml {z}) {m^ (z))*] 
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where ojm = We then apply the Fourier transform to move into the frequency do¬ 
main, where dtm^ z) ^ {Hq, I; z), such that the energy relaxation rate (1/Ti)* 

for magnetization oscillations along the x-axis is 


where 


/ 1 V _ 


(23) 


A',. := 


oo S-\-d _ 

J dx J dz {m^{z)) {x, z) 

-OO d 


S+d 


(24) 


J dx J djlm|m;(j)(m;(z))' 


is a dimensionless inductive coupling parameter. In the limiting case of the n = 0 (i.e., 
uniform) mode with a uniform excitation held due to current howing only through the 
waveguide center conductor, and an inhnitesimal spacing between the waveguide and the 


sample, we have Kq = 6/Aw. Substituting Eq. (15) into Eq. (24), we obtain the general 
result 


Kn 


OO 6^d _ 

f dx f dzqn {z) {x, z) 

— OO d 

5-\-d 

e J dz\qn{z)f 

d 


(25) 


with e = |m^| / \ fh^\. 

Since the energy dissipation rate for the case of radiative damping is anisotropic, it must be 
generally treated in the damping tensor formalism, where the Gilbert damping torque T is 
given by 


Tk = eijkaijTTiiidtm).. (26) 

The equation of motion is 

dtifi = —jfio'di X -I- T (27) 

and m = M/Mg is the normalized magnetization. For the coordinates in Fig the 
only nonzero radiative damping tensor components are a^x and ayx- For the perpendicular 
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FMR geometry, the relationship between the energy relaxation rate and the Gilbert damping 
components is 


Ti 


— OLy^cc.^:r. 


(28) 


and 


Ti 




(29) 


where and Uy are the respective stiffness frequencies, defined as 


OJi : = 


7 




Ms dmi 

and Ujn is the magnetic free energy function. The frequency-swept linewidth Ao; = 
7 /ioAif, where Aif is the field-swept linewidth in Eq. 10. is given by 


Aw = 




^zx^x ^zy^y 


(31) 

(32) 


For perpendicular FMR, = Uy = w, and the specific case of anisotropic radiative 
damping, a^x = <^zy = 0, and we obtain 


= 




K„ 


(33) 


2w \Ti / ^ 2Zo 

and This is in contrast to the case of isotropic damping processes, such 

as eddy currents and intrinsic damping, where we obtain Aw)^° = 2a'^°uj instead. Thus, 
the net damping due to the sum of anisotropic radiative damping, and any other isotropic 
processes, is given by 

^rad 

(34) 


„ „ int I „ eddy , 

Oin — ® "I-7 


where is the damping parameter in Eq. (11) for the field-swept linewidth. 


We use a vector network analyzer (VNA) to measure the two-port R-parameter matrix 
element for the nth spin wave mode, AS*^^. The matrix element is defined as the ratio of the 
voltage induced in the waveguide by the nth spin wave mode 14i(-ffo) in an applied magnetic 
field Hq, and the excitation voltage Gm 
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:= 


VnjHo) 
K;n ' 


( 35 ) 


If we model the reactance of the nth spin wave mode as nothing more than a pnrely indnc- 
tive element of indnctance in series with an impedance matched transmission line, and 
if we assnme the sample indnctance is mnch smaller than the transmission line impedance, 
we can approximate as 


AS^\Ho) = 


iuLniHo) 


(36) 


where L„(hro) = {Ho, I) /I. 


We dehne a normalized, held-independent mode-indnctance as 


Lr, . 


Ln{Ho 


xn{Hoy 

and a dimensionless, field-independent mode-amplitnde A^, 


An, : — 


iuLn 


snch that 


(37) 


(38) 


ASl\H„) = -A„x„(Ho). (39) 

Thns, An is the dimensionless amplitnde parameter that we obtain when htting data for 


ASy{Ho). By use of Eqs. 18, 38, and 37, we can rewrite the mode-amplitude as 


An • ‘I 


. Ufiol 


fS+d N 2 

dx I dzqn{z)h^{x, z) 


(40) 


WiZ„ v-oo 

Remembering that the normalized mode inductance has a factor identical to the numer¬ 


ator of Eq. (25), we can rewrite the radiative damping in terms of the normalized mode 
inductance. 





where 


Vn ■= 


5 

dz\qn{z)\^ 


(41) 


(42) 
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We emphasize that Eq. (41) is a very general result, regardless of the details of the 
excitation held prohle . Thus, even if the held prohle is highly non-uniform due to the 
combination of eddy current and capacitive coupling ehectd^^®, there should still be a hxed 
scaling between the radiative damping and the normalized inductance. 


In the case of the uniform mode, t] = 1/A and Eq. (41) reduces to 


a 


rad 


Iz„ 


(43) 


However, for a sinusoidal mode of the form 


qn{z) = cos 


(2n -|- 1) Tiz 
25 


(44) 


that is expected in the case of a pinned boundary condition at one interface and an open 
boundary condition at the other interface, as shown in Fig. (b), we obtain rj = 1/2 and 


a 


rad 


ojm 


I - 2Z„' 

±jn ^ 

For the case of the wavenumber values extracted from the data shown in Fig. (c) for a 
200 nm Py him, we can determine the value for and the degree to which it can vary with 
mode number. We use the following form for the spin wave prohle: 


qn{z) = cos {knz) (46) 

consistent with our assumption, when extracting kn from our PSSW data, that an un¬ 
pinned boundary condition applies to only one of the interfaces, i.e. at z = 0. Using these 
extracted values for the wavenumber, we obtain values for rjn shown in Fig. 

We see in Fig. that the variation in with varying mode number is less than 10%. 
Thus, to within hrst order, we can treat as a constant for the purposes of htting our data, 
i.e. rjn = r]. 


B. Derivation of 

For the derivation of the eddy current damping Oeddy uniform magnetization dynamics 
are assumed. The notation stays the same as for the radiative damping. 

Then the total hux passing through the magnetic him is 
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Figure 9; Dependence of calculated values for rjn on mode number for the case of the spectral data 
presented in Fig. (c) and[^ (a). Wavenumbers are extracted from those data via the procedure 
outlined in the main text of the paper, based upon a model with a single surface with interfacial 
anisotropy, and an unpinned boundary condition at the other interface. Within the context of that 
particular model, and the expected quadratic dependence of spin wave resonance frequency with 
wavenumber that it produces, we see that has a weak dependence on mode number, justifying 
our presumption that r]n can be treated as a constant for the purposes of fitting the damping data. 


dt^ = fio^6{dtrn)^, 

where = x ■ dtfh. The electrical power dissipated by the eddy-currents is 


(47) 


^ 1 

2 f 2pi A 

ys.aw J 


--:- \{dtm) 


(48) 

(49) 


with ^eff := C6/2, where 0 < C < 1 is a phenomenological parameter that accounts for 
details of the non-uniform eddy-current distribution in the ferromagnet. Analogous to the 
derivation of the radiative damping, we now need the energy of the magnetic excitations. 
The number of magnons in the system is given by 


N = 

’mag 








( 50 ) 
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Thus, the total magnon energy is 


Ejnag i.6 

= —\dtmfWi6. 

-yu 

The rate of energy dissipation is then given by 

^ _ ^ind _ C 'yufilMsS'^ 

Ti -Emag 8 p \dtrh\^ 

The maximum energy decay rate occurs when {dtrh)^ = \dtm\, in which case 


(51) 

(52) 


(53) 


where the superscript indicates that this is the maximum decay rate for magnetization 
oscillations along the x-axis. For the case of a perpendicular applied held sufficient to 
saturate the static magnetization out of the him plane, the damping process is isotropic, 
i.e.. 


1 y _ 

Ti) 8 p 


Therefore, analogous to Eq. (30), the frequency-swept linewidth Aon is simply 


(55) 





then 


16 p 


(56) 
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